Swimming speeds of filaments in nonlinearly viscoelastic fluids 

Henry C. Fu/'Q Charles W.Wolgemuth,^'Q and Thomas R. Powers^ 

^Division of Engineering, Brown University, Providence, RI 02912 
^Department of Cell Biology, University of 
Connecticut Health Center, Farmington, CT 06030 
(Dated: June 12, 2008) 

Abstract 

Many microorganisms swim through gels and non-Newtonian fluids in their natural environ- 
ments. In this paper, we focus on microorganisms which use flagella for propulsion. We address 
how swimming velocities are affected in nonlinearly viscoelastic fluids by examining the problem of 
an infinitely long cylinder with arbitrary beating motion in the Oldroyd-B fluid. We solve for the 
swimming velocity in the limit in which deflections of the cylinder from its straight configuration 
are small relative to the radius of the cylinder and the wavelength of the deflections; furthermore, 
the radius of the cylinder is small compared to the wavelength of deflections. We find that swim- 
ming velocities are diminished by nonlinear viscoelastic effects. We apply these results to examine 
what types of swimming motions can produce net translation in a nonlinear fluid, comparing to the 
Newtonian case, for which Purcell's "scallop" theorem describes how time-reversibility constrains 
which swimming motions are effective. We flnd that the leading order violation of the scallop the- 
orem occurs for reciprocal motions in which the backward and forward strokes occur at different 
rates. 
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I. INTRODUCTION 



Bacteria and other microscopic swimmers live in a world of low Reynolds number, in 
which viscous effects dominate over inertia [l, 2|. In this regime, the drag forces acting on a 
swimmer are much greater than the forces necessary to achieve the observed accelerations of 
these swimmers. At the limit of zero Reynolds number, the sum of the net drag force and any 
other external forces, such as buoyancy, must vanish. Therefore, the physics underlying the 
swimming strategies of microorganisms is markedly different from that of large organisms 
such as fish. For example, in the absence of inertia, the diffusion of viscous shear waves is 
essentially instantaneous. Thus, the flows induced by a swimming creature's motions start 
and stop immediately when the creature's movements start and stop. This phenomenon 
is closely related to the "scallop theorem," of Purcell j2|, which states that time reversible 
strokes in a Newtonian fluid, such as water, cannot produce net swimming because there 
is an equal and opposite net translation for corresponding forward and backward stroke 
motions. As a result, swimming microorganisms use strokes with a clear time direction, 
such as the traveling wave patterns on the flexible flagella of mammalian sperm, or the 
rigid-body rotation of the helical flagellar filaments of Escherichia coli. 

Early theoretical work elucidated the basic principles of swimming at zero Reynolds 
numbers in a Newtonian fluid. For example, Taylor calculated the swimming speed of 
an infinite two-dimensional sheet as a function of the amplitude, wave number, and wave 
speed of a traveling wave on the sheet f3|]. He also did the same calculation for an infinite 
cylindrical filament jj] . These calculations showed that the swimming speed is proportional 
to the wave speed, as expected from dimensional analysis. Since the shape and speed of 
the waves on the sheet or filament in these calculations are independent of the load, the 
swimming speed is independent of the viscosity of the fluid. In other work, various groups 
have calculated how viscous forces and forces internal to a flagellum determine the shape of 

SkZI or actuated by internal motors distributed along 



flexible filaments actuated at one end 
the filament js-ll]. 



All the calculations just mentioned apply to Newtonian fluids. However, swimming mi- 



croorganisms often encounter elastic fluids. For example, mamma! 



through viscoelastic mucus in the female reproductive tract 
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cobacter pylori must swim through mucus lining the stomach 15|; and spirochetes must 



lan sperm must swim 
J]; the bacterium Heli- 
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burrow through connective tissue in the course of infection [l^. The theory of swimming in 
non- Newtonian fluids is much less developed than the Newtonian case. Fulford, Katz, and 

n 

Powell extended the resistive force theory of Gray and Hancock [17] to the case of a linear 
Maxwell model for viscoelasticity, and found no change in the swimming speed compared to 
the Newtonian case More recently, Lauga emphasized that since swimming speed de- 
pends nonlinearly on wave amplitude, a realistic treatment of swimming in a non- Newtonian 
medium should include a nonlinear model for the viscoelastic medium fio!]. Lauga inves- 
tigated several nonlinear fading memory models and found that for a given wave pattern 
on a two-dimensional sheet, viscoelastic effects slow the swimmer relative to its speed in a 
Newtonian liquid jiol - 

In our previous work, we calculated the swimming velocity of an infinite filament with 
a prescribed waveform in the Upper Convected Maxwell Model 20|, a simple nonlinear 



model for a viscoelastic material with fading memory 2l|. We also studied the effect of 
the viscoelastic medium on the shape of a beating filament described by a simple sliding 
filament model js], and showed how altered beating shapes affect swimming velocities. In 
the current paper, we give a more detailed account of our calculation of the swimming 
speed of an infinite filament with a prescribed waveform. We refine our previous calculation 
by including the effects of solvent viscosity with the Oldroyd-B model for a nonlinearly 
viscoelastic polymer solution. As an application of our calculations, we examine how the 
scallop theorem is spoiled in a nonlinearly viscoelastic fluid at zero Reynolds number. In 
this paper we consider prescribed shapes of beating patterns. In actual swimmers with 



beating filaments, the shape is determined by the interaction o: 



hydrodynamic forces, which we address in another publication Here we take the point 



filament structure with 



of view that the beating pattern is already known (for example, from direct observation) 
and ask what the consequences are for swimming properties. We restrict our discussion to 
time-periodic motions and associated flows, corresponding to repetitive swimming actions 
that can be employed to move an arbitrarily large distance by a self-contained organism. For 
these time-periodic motions, the associated flows will have both time-periodic and constant 
components. Only the constant, time- averaged component produces net translation over 
many cycles. 

We close this introduction with a summary of the structure of this paper. In Sec. HTl we 
review the Oldroyd-B model and how it naturally arises for a material with fading memory. 



3 



Then, in the next section, we calculate the flow fields induced by a helical traveling wave 
on an infinite cylinder. We work to second order in the defiection of the cylinder, and 
deduce the time-averaged swimming speed. In Sec. IIVI we explore the ramifications of our 
calculation for the scallop theorem, and show quantitatively how a reversible stroke pattern 
leads to propulsion. Section |V] is the conclusion. Finally, in the appendix we calculate 
the viscoelastic force per unit length to first order acting on a filament with transverse or 
longitudinal (compressional) traveling waves. 



II. NONLINEARLY VISCOELASTIC FLUIDS 



Since we consider the zero Reynolds number regime, inertia is irrelevant and the dynamics 
of an incompressible medium is governed by force balance, — Vp + V • r = 0, where p is 
the pressure associated with incompressibility, V ■ v = 0; v is the fiow velocity; and r is 
the deviatoric stress tensor. The particular character of a medium is determined by the 
constitutive relation, which relates the stress to the strain and rates-of-strain. In this paper 
we consider dilute polymer solutions, in which the stress relaxes to the stress of a Newtonian 
viscous fiuid over a time scale A. A simple constitutive relation for a dilute polymer solution 
is the single-relaxation-time Maxwell model: 

nj + A-^ = 7]^ij, (1) 

where i and j label Cartesian coordinates, rj is viscosity, and = dvi/dxj + dvj/dxi is the 
strain rate. This constitutive relation is valid to first order in deformation. However, the 
swimming velocity of a prescribed motion is second order in the deformation 3|, To see 
why, consider a traveling wave of amplitude h on an infinite cylindrical filament (Fig. [1]). 
The swimming velocity may be represented as a power series in h. Replacing h with —h 
changes the sign of the first order contribution to the velocity. On the other hand, changing 
the sign of h amounts to translating the filament by half of a wavelength, and therefore yields 
no change in the time-averaged swimming velocity. Thus we conclude that the swimming 
velocity is quadratic in h to leading order, and therefore nonlinear corrections to Eq. ([T]) 
will be important for the study of swimming in a viscoelastic medium jl^. 

It is useful to recall why the constitutive relation ([1]) is inapplicable for large deformations. 
First, note that the Maxwell constitutive relation does not respect Galilean invariance. This 



4 



defect is easily fixed by replacing the partial time derivative of the stress, drij/dt, with the 
material derivative, drtj/dt + v • Vrjj. However, this modified constitutive relation is still 
unsuitable for working beyond first order in deformation because it is sensitive to the overall 
orientation of the material. For example, consider a steady-state flow, for which ([T]) predicts 
that Tij = rj'jij. Now place the medium on a slowly rotating turntable [20i]. Since viscous 
forces are much greater than the fictitious forces arising from rotation, the relation between 
stress and strain rate in the rotating material is the same as that in the stationary material, 
to an excellent approximation [22i]: t-j = 777^'^, where the prime denotes the quantities for 
the rotating material. On the other hand, the components r/^- are time-dependent, because 
= Rii^RjiTki, where Rij is a time-dependent rotation matrix. Therefore, Eq. ([T]) cannot 
hold for the primed quantities, even with the partial time derivative replaced with the 
material derivative. 

In general, changes in the components t^j due to motion of the basis vectors should not 
enter the constitutive relation. The formulation of rheological equations of state requires a 
time derivative that does not depend on the choice of basis [23] . Any constitutive relation 
that reflects the small magnitudes of fictitious forces must be formulated in terms of time 
derivatives of stress or strain that transform homogeneously under change of coordinates. 
There are many different ways to construct a time derivative of the stress tensor that is 
insensitive to the rate of change of the basis vectors. In the following we present one way to 
do so which leads to the nonlinear viscoelastic model used in the rest of the paper. 

Due to the small magnitudes of fictitious forces, the constitutive relation should be insen- 
sitive not only to global rotations such as that of the turntable described above, but also to 
local rotations and translations of fluid elements. One way to enforce this more general con- 
dition is to express the time derivatives in terms of coordinates which track fluid elements. 
Suppose {^^,^^,^^} are convective coordinates, with r(^^, t) the position at time t of 
the particle with Cartesian coordinates {^^, C,^} at the initial time. Convective coordinates 
are also known as Lagrangian coordinates. The basis vectors associated with convective co- 
ordinates are gj = dr/dC,^, and the stress takes the form r = r^^gi ® gj, or more simply 
T = gjT^-'gj. For Cartesian coordinates, whether an index is raised or lowered has no sig- 
nificance; however, for general coordinates, raised indices signify contravariant coordinates, 
and lowered indices signify covariant coordinates [20]. Since dgi/dt = dv/d^^ = (gj ■ V)v, 
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we have 

The second two terms in Eq. (|2]) arise from the the dependence of the basis vectors on 
time and therefore they cannot enter the constitutive relation. To ehminate the spurious 
time-dependence, define the upper-convected time derivative using convective coordinates: 



" dt 

In coordinate-free form, the upper-convected derivative is 



g.- (3) 



T= + V ■ Vr - T ■ Vv - (Vv)^ ■ r. (4) 
at 

It follows directly from Eq. ([3]) that the components of r in a general coordinate system 
transform homogeneously even under a time-dependent change of coordinates. 

The upper-convected time derivative is just one way of many to construct a time derivative 
of the stress tensor that is insensitive to the rate of change of the basis vectors. For example, 
if we consider the covariant components Tij of the stress tensor in the convected-coordinate 



basis, a different natural derivative emerges: the lower-convected derivative [2J]. The choice 
of one particular convective derivative over another cannot be justified on symmetry grounds 
alone, but must be made using information from experiments or the continuum limit of a 
microscopic model. 

For example, a commonly used model for dilute polymer solutions is the Oldroyd-B 
model 20|]. We will use this model in this paper. In the Oldroyd-B model, the polymer 
component is represented by the Upper Convected Maxwell model, which is the Maxwell 
constitutive equation ([1]) with the time derivative replaced by the upper-convected time 
derivative. The appearance of the upper-convected time derivative, rather than some other 
form insensitive to rates of change of basis vectors, can be justified by appealing to a mi- 



croscopic model of polymer stresses arising from Hookean dumbbells [25|]. The solvent is 
represented by the usual Newtonian constitutive relation for a viscous fluid. Thus, the 
deviatoric stress for the solution is 

r + Xt= 7]j + Ar^s 7, (5) 

where the total viscosity rj = rjp + r]s is the sum of polymer and solvent viscosities. Note 
that if a patch of fluid undergoes a spatially uniform rotation or translation, due to the 



covariant nature of this equation, within the patch the constitutive relation is insensitive to 
the motion. In (|5]), if the solvent viscosity is ignored (r/s = ), the Oldroyd-B fluid reduces to 
the Upper Convected Maxwell fluid, and if there is no polymer [f] = f]s) the Oldroyd-B fluid 
reduces to a Newtonian fluid. The model for the Oldroyd-B fluid predicts realistic elastic 
and first normal stress effects as long as elongational flows are not large; however, it does 
not capture shear-thinning or yield-stress behaviors observed in some non-Newtonian fluids. 



The Oldroyd-B model has been used to describe Boger fluids [26|. We use the Oldroyd-B 
fluid as an example. The choice of constitutive relation ultimately depends on accurately 
modeling the properties of particular real-life fluids and flows. We note that in the case of 
an infinite sheet, Lauga examined a number of nonlinear models, including Finite Extensible 
Nonlinear Elastic models which remedy the limitations of Oldroyd-B in elongational flows. 



and found that the swimming speed was not affected by the choice of constitutive model [1 

III. SWIMMING VELOCITY OF FILAMENTS IN VISCOELASTIC FLUIDS 

In this section we describe the solution to the problem of an infinitely long cylinder 
with prescribed beating motion in a viscoelastic fluid. We obtain the swimming velocity 
in a nonlinearly viscoelastic fluid, and find that in general it is decreased relative to the 
swimming velocity in a Newtonian fluid. Our results are different from the results of Fulford, 



Katz, and Powell 18|, who used the constitutive relation ([T]) to find that a filament with 
prescribed beating pattern swims at the same speed in viscous and linearly viscoelastic fluids. 
Our results are in broad agreement with those of Lauga, [L9\ who solved the corresponding 
problem of an infinite sheet with traveling-wave displacements. 

Figure ([1]) shows a snapshot of the swimming filament. The beating motion is represented 
by a transverse traveling wave. This motion induces flow. To calculate the swimming 
velocity, we work in the frame of the swimmer and calculate the z-component of the flow far 
from the swimmer. As argued above, the swimming velocity is second order in the amplitude 
of the deflections of the filament; therefore we must calculate the ^-component of the flow to 
second order. We calculate the average speed over one period of the motion of the filament. 
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FIG. 1: A section of an infinite filament with prescribed beating pattern h{s,t). 

A. First order solutions for a cylinder with prescribed beating pattern 

In this subsection we define the deformation of the filament, prescribe the boundary 
conditions, and solve for the flow velocity to first order in the deformation of the filament. 
We consider a cylinder with radius a, with its axis lying along the z direction. For the 
undisturbed, straight cylinder, material points on the surface are parameterized by s, the 
distance along the axis, and the angle ip (see Fig. [1]). Note that these coordinates are 
convective coordinates for the filament surface. For a transverse planar traveling wave, the 
deformation of the cylinder is given by prescribing the position sz + h{s, t)x of the centerline 
of the cylinder. Material points on the surface of the cylinder are parameterized by 

r{ip,s,t) = ^[h{s,t) + a cos ip] + y a sin ip + zs; 

h{s,t) = Re{h exp{iqs — iut)} . (6) 

Because the displacement is solely in the x direction, the filament must be extensible. Enforc- 
ing inextensibility leads to motion of the filament surface in the z direction with amplitude 
of order h"^; this type of longitudinal motion would give rise to corrections to the swimming 
velocity of order h'^, higher order than the results we report in this paper. For technical 
ease, we solve for arbitrary motion by decomposing into a linear superposition of traveling 
waves; therefore it is sufficient to assume that the prescribed shape change in Eq. (|6]), is 
given by a single propagating mode. As we will show, the swimming velocity, although 
second order in h, has a form amenable to this decomposition into traveling waves. We 
also introduce notation for variables with tildes that we use in the remainder of this paper: 
A = Re{Aexp {iqs — iut)}. 

Assume that the displacement of the filament from its undisturbed state is small compared 
to the smallest wavelength of the cylinder motion h{s, t) and the cylinder radius a. This 
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assumption allows us to take advantage of the cylindrical symmetry. Later, we will find it 
convenient to further demand that ga ^ 1. Using cylindrical coordinates, with basis vectors 
z and 



f = X cos ip + y sin %l) (7) 
(f) = — xsin-?/; + y cos-?/), (8) 

denote the position of points on the surface by (rb,0b,^b)- (The subscript b stands for 
"boundary.") To first order in h (see Fig. 

Tb = a + h cos ip 

0b = 'ip sin ip 

a 

Zb = s. (9) 
For no-slip boundary conditions on the filament surface, 

(10) 



dr 

v(r(V^,s,t)) = — 



■4>,s 

Since our goal is to calculate the swimming velocity to leading order, we must expand the 
boundary condition (|TOl) to second order in the dimensionless displacement qh. To all orders 
in qh, the right-hand side of Eq. (fTOj) equals /ix, where here and elsewhere the dot signifies 
the partial time derivative. Since the leading order contribution to the fiow velocity is first 
order in h, we need only expand the left-hand side of Eq. ( [T0|) to first order in the argument 
of the fiow velocity: 

W + h^-'^<,^,s,t = h^- (11) 
The fiow velocity at a fixed point also has an expansion in powers of qh: 

v = vW + v(2) + .... (12) 

Thus the boundary condition to first order is 

v^^\a,ip,s,t) = h± (13) 



To obtain the forces and swimming velocity, we need to solve the Stokes equations men- 
tioned above, V ■ r — Vp = and V ■ v = 0, where r is now the deviatoric stress tensor 
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FIG. 2: A cross-section of the infinite filament, showing the undeformed (dashed hne) and current 
configuration (solid line), and the relation between (p]^ and the convective coordinate ■0- 



that satisfies Eq. To simplify the discussion of the equations of motion for the fluid, we 
nondimensionalize by measuring length in units of 1/q, time in units of 1/uj, and pressure 
and stress in units of rju. Hereafter in this and the next section, we use the same variables 
to denote the nondimensional quantities, unless explicitly stated otherwise. At linear order, 
all time- dependent quantities are proportional to exp(— it).To first order, Eq. (|5]) is 

l-iDe2r(i) 

where the Deborah numbers De = cuA and De2 = Deris/rj. The first-order flow obeys 
V ■ v(i) = and 

y2 _ IziiHiv (p(i)e-^(*-^)) = 0. (15) 

^ ^ 1 — iDe2 

The flrst-order boundary condition f|T3|) reduces to 'v^^\a,'tp) = — i/ix, where v^^^^ = 
Re{v^^^(r, (f)) exp(i2; — it)} in accord with the convention of Eq. ([6]). 



The solutions to these equations in cylindrical coordinates are known J, l27|]. Note that 



althougli there is no solution to the Stokes equations for dragging a rigid inflnite rod through 
fluid [28| , the flnite wavelength together with the rod radius provide a cutoff in our problem. 
The flrst order results for the pressure and velocity fleld are presented in Table [B In the 
table, KqIt) and Ki{r) are modifled Bessel functions of the second kind. These results are 
valid up to leading order in an expansion in 1/ log(a). The flow flelds are the same as those 
in the Newtonian case ^|, but the pressure has an additional complex factor. 

These flrst order solutions are necessary for the calculation of the swimming velocity of 
the fllament, which we present in the next subsection. These solutions are also sufficient for 
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p 

Vz 

aA 
B 
C 



— i/icos (j)AKi{r) 
-\hcos(l)[aArKi{r) + BK2{r) + CKQ{r)] 

-ihsm(P[BK2(r) - CKo{r)] 
-h cos (/) [aArKo{r) + {B + C - aA) Ki{r)] 



^Ko{a) + aKi{a) 



1 , Koja) Kojaf 



+ 



27^2 (a) i^i(a)2 



r 



-aAaKi{a) / [2K2{a)] 
[1 - aKi{a)aA/2] /Ko{a) 



TABLE I: Pressure and velocity field for a cylinder with a transverse traveling wave; a = (1 — 
iDe)/(l -iDea). 

calculating the lowest order effects of viscoelasticity on hydro dynamic forces and the shapes 



of flexible swimmers 



In appendix 1^1 we calculate the force per unit length acting on 



filaments with transverse and longitudinal traveling waves. These results are consistent with 
the resistive force theory of Fulford, Katz, and Powell 18 1. 



B. Velocity of a cylinder with prescribed beating pattern 

Turning now to the second-order calculation, we expand the boundary condition (flOl) to 
second order: 

To isolate the nonlinear terms of the equation of motion, it is convenient to rewrite the 
convected time derivatives as 



T = dtT + T 

i = da + G. 



(17) 
(18) 



The second-order equations are 



:i + De5i)r(2) = (1 + De2dt)j^^^ - DeT^^) + DeaG 



(2) 



(19) 
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or, invoking force balance — Vp^^-* + V ■ r^^^ = 0, 

(1 + Dedt)Vp^^^ = (1 + De29t)VV(2) - DeV ■ T^^^ + DeaV ■ G^^^. (20) 

Note that T*^^) and G^^-* only have contributions from first order flows and stresses. 

We must solve Eq. ( I20l) for Vz ■ The time dependence of Vz is most easily expressed by 



expanding in modes e*™(^~*). In our first order calculations, only m = 1 appeared. Second 
order involves m = 0, ±2. Similarly, the azimuthal dependence of velocity and stress fields 
can be expanded in modes of cos(n0) and sin(n0), with n = and n = ±2 appearing at 
second order. To find the swimming velocity, it is only necessary to solve for the time- 
and 0-averaged ^-component of the velocity field, which corresponds to the net flow in the 
2;-direction far from the filament. This net flow is the m = 0, n = mode. In the m = n = 
part of the second-order equation ( |20i) . all terms with (9^, dt, and dz drop out, leaving 



Tn(2)(m=n=0) _ r»(2)(m=n=0) 



where Trz can be expressed in cylindrical coordinates as 20 1 



(21) 



Trz = VrdrTzr+VzdzTzr + — d^Tzr-TrrdrVz-Tzr{dzVz + drVr)-TzzdzVr- — d^Vz-—d^Vr. (22) 

ly ly ry 

A similar expression holds for Grz-, with the components of r replaced by the corresponding 
components of 7. Examination of the real part of Grz reveals that every term is proportional 
to an odd power of sm{z — t) or cos(2;— t); therefore, the average of Grz over a period vanishes. 
Integrating Eq. (I2T]) over r yields 

^(2)(™=„=o) = Dey drT42){™="=0) - f/, (23) 

where t/ is a constant. The complete expression for Tr^^ in terms of a and the modified Bessel 
functions -ft'o(^) and Ki{r) may be constructed from the entries of Table [H To simplify the 
resulting expression, it is convenient to expand the coefficients aA^ B, and G in powers of 
a. In this expansion, we only keep the leading logarithms. For example, we approximate 
Ko(a) ~ log2/a — 7 ~ — log a, where 7 is the Euler constant. Since we use the boundary 
condition for vi^^ at r = a to determine the swimming velocity, we next expand 7"J^)(m,n-o) 
in powers of a to find 

y(„o) ^ logr 

2 l + De^ rlog'a ^ ^ 
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The second order velocity in the m = 0, n = mode is thus 



|/iplogVDe - De2 
2 log^ a 1 + De2 

To determine the constant U, average the boundary condition (fT6|) over t and to find 



^ ^ 21og% 1 + De2 ^ ^ 



^(2)(m=n=0) I 



-/i COS (pdrV^p H — sin (pd^v^^^ 



{m=n=0) 



hY/2. (26) 



The swimming velocity is therefore 



|/ip l + DeVr/ 

(7 = — ^ — , (27) 

2 l + De^ ^ ^ 



where we have reinstated the dimensions and used De2 = Deris/i]. The filament swims in 
the same direction as it would in a Newtonian fiuid, opposite the direction of the traveling 
wave, but with a reduced speed, since rjs/r] < 1. This result is valid up to leading order in 
1/ log a, and to second order in h. 

In the appendix we describe the first-order solutions for a longitudinal wave. The same 
steps just outlined lead to a swimming velocity of 

\h\ H + De\/rj 
Uiong = ujq— Y+B^' ^ ' 

For a longitudinal wave, the filament again swims in the same direction as it would in a 
Newtonian fiuid, along the direction of the traveling wave, and with a reduced speed. To see 



why an object with 
wave, see reference 



ongitudinal propagating waves swims in the direction of the propagating 



29| 



To apply these results to filaments with generic beating patterns, we must consider the 
swimming velocity of a cylinder with a beating pattern described by a superposition of 
transverse or longitudinal traveling waves, 

h{s,t) = J2K,^e'^''"^'\ (29) 

written in dimensional form. Since the first-order equations are linear, the first order fiow is 
a superposition of fiows corresponding to the individual Fourier components hq^^j- As above, 
the swimming velocity is affected only by TrT~^~^^ and G'^~^~^\ Due to the averaging 
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over t involved in obtaining the m = mode of Trz, the first order flows with different 
wavelengths do not interact, yielding 

rJr="=°^ = E ^r^r="=°^ ^) , (30) 

q,u> 

where Tr^~'^~^\q,uj) is TrT~"'~^^ calculated for a single traveling wave. This separation of 
modes holds for the swimming velocity, and we obtain the dimensional result 

u = -l^^^j^i^\KJ . (31) 

Note that the factor of two difference between Eqs. fl271) and fl3ip arises since the single wave 
cos(s — t) has /i = 1 in Eq. (Q, and hu = /i-i-i = 1/2 in Eq. fl29|) . Finally, we note that 
for three-dimensional transverse motion, 

h{s, t) = J2 [e.hlJ^'^--'^ + eyhU^'''-'^'^] , (32) 

q,U} 

the two polarizations decouple upon averaging over to obtain the n = mode, leading to 

^ = - i + iooxy U^'-' + ' J • ^^^^ 

IV. RAMIFICATIONS FOR THE SCALLOP THEOREM 

As mentioned in the introduction, the kinematic reversibility of Stokes equations leads 
to the scallop theorem for viscous fluids, which states that reciprocal motion of a body 
leads to no net translation. A reciprocal motion is one which is time reversible, up to a 
reparameterization in time. In this section, we use the swimming velocity in an Oldroyd 
B fluid, Eq. flHTl) . to examine how the scallop theorem breaks down in viscoelastic fluids. 
Equation fl3T|) is only valid to second order in displacements, and all of our statements 
have the same limitation. Due to this limitation, motions which we identify as having 
zero swimming speed may have nonzero swimming speeds at higher order; nonetheless, the 
reciprocal motions that have nonzero swimming speed in Eq. fl3ip serve to identify the 
leading order violations of the scallop theorem. 

We start with some specific classes of beating motions that satisfy the hypotheses of the 
scallop theorem. In these examples the motion is confined to a plane, but the conclusions 
can be generalized to three dimensional motion. Consider a strictly time-reversal invariant 
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beating motion, h{s,t) = h{s,—t), where h is a transverse deflection, and, without loss 
of generality, the origin has been chosen as the midpoint of the period. In this case, the 
Fourier coefficients satisfy hg^^^ = hq_^, and therefore the sum in Eq. (ETl) vanishes. Note 
that the same argument can be made for space reflection symmetric beating patterns as well. 
Next, consider a beating motion which has effectively only one degree of freedom: h{s,t) = 
f(t)g{s). In this case, the Fourier coefficients are also factorizable, so hg^^ = f{Lo)g{q). Since 
both f{t) and g{s) are real functions, f{uj) = f*{—uj) and g{q) = g*{—q), and therefore 
m,ajp = As a result, the summand in Eq. (!3T]) is odd in co, and the swimming 

speed vanishes. More generally, any motion with = \hq-uj\, which corresponds to a 

superposition of standing waves, leads to no net translation. 

These examples demonstrate how time-reversible motions lead to zero swimming speed 
at leading order in an Oldroyd B fluid. In a Newtonian fluid at zero Reynolds number, more 
general reciprocal motions, such as a fast power stroke followed by a slow return stroke, also 
lead to zero swimming speed. To see how this conclusion does not apply to the Oldroyd B 
fluid, it is convenient to write the swimming speed U in the time domain and in real space. 



Any reciprocal motion h{s, t) with period T may mapped to a strictly time-reversal invariant 
motion by a reparameterization ti = F{t), where F is monotonic, F{0) = 0, and F{T) = T. 
For convenience we have chosen the origin as the beginning of the period. In the Newtonian 
case, A = 0, and the kernel in Eq. becomes a Dirac delta function, M{t,t') = 6{t~^ —t'), 
with the interpretation that the singularity of the delta function lies within the region of 
integration. Inspection of the resulting expression for U reveals that it is reparameterization 
invariant, and we conclude that our formula for the swimming speed obeys the scallop 
theorem for general reciprocal motions in a Newtonian liquid at zero Reynolds number. 
However, when A 7^ 0, the kernel M{t,t') spoils reparameterization invariance, and the 
speed need not vanish. The physical interpretation is that the distance covered per period 
in a non-Newtonian fluid depends on the rate of motion as well as the sequence of shapes. 
For example, consider a helical wave which propagates forward by one wavelength in the 




(34) 



where the memory kernel M(t, t') is the operator 




(35) 
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first third of tlie time period, and tlien backwards by one wavelengtli in tlie final two-thirds 
of the period: 

{bcos{s — 3t)ex + bsin{s — 3t)ey < t < 27r/3 

(36) 
-6cos(s + 3t/2)e^ - 6sin(s + 3t/2)ey 2n/3 < t < 2n 

For this motion, the Fourier amplitudes for u = 3 are /if 3 = h^^ _^ = ih\^ = —ih^_i = b/6. 
For w ^ 3, we have h'^^ _^ = {hf^^)*, = -i/if_^, and hl^ _^ = {hl^)* , where 

_ 96 exp (-i7r(x;/3) sin {7iu/3) 
'''^ 47r (a;-3)(a; + 3/2) ' 

All other hg^^^ are zero. Note that these coefficients do not obey the criterion \hq^i^\ = 

\hq-uj\. Inserting these coefficients into Eq. (!33|) and using contour integration to evaluate 

the resulting sums 



U 



30| yields the dimensionless speed 

Vr(4 + 45De2 + 81De^) + 18De(2+9De^)sinh|5, 



27De%^ 



1+2 cosh ^ 



(38) 



r] J n (4 + 45De^ + SlDe^)^ 

where De = 271X1 T and T is the dimensional period of the motion. 

The swimming speed vanishes in the Newtonian limits of no polymer, rjs = rj, or zero 

relaxation time, De = 0. For these cases, the swimmer moves in the negative direction 

during the first third of the period, and in the positive direction by the exact same amount 

during the last two thirds of the period. In the viscoelastic case, the swimmer moves more 

slowly during the first third of the period when the wave frequency is high. Thus, the 

scallop theorem is rendered inapplicable and the net motion of the swimmer is in the positive 

(forward) direction. The formula (138|) for the swimming speed is complicated due to memory 

effects: the speed of the swimmer during each stroke depends on its motion during previous 

strokes. However, we can make a simple estimate for the swimming speed by disregarding 

these memory effects and using Eq. flTTj) to estimate the speed during each stroke of the 

waveform of Eq. : 

_ [ 1 + (9/4)De^(r/,/r^) 1 + 9De^(r^,/r/) ] 

[ 1 + (9/4)De2 1 + 9De' \ ' ^^^^ 

Figure [3] shows the exact second-order ( l38i) and estimated ( l39l) swimming speeds for the 

waveform (1361) . with dimensional wavelength 27i/q. 

Helical waves provide a convenient experimental system to explore the fate of the scallop 

theorem in viscoelastic fluids, since the backward and forward travehng waves could be easily 
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FIG. 3: (Color online) The swimming velocity of the reciprocal motion of Eq. (j36p versus De = 
27rA/T, for rjs/rj = 0.75. The (red) solid curve is the exact second-order result of Eq. (j38|) : the 
(blue) dashed curve corresponds to the estimate of Eq. 



actuated by rotating a rigid helical filament wrapped in a flexible rubber sleeve |^]. For a 
cylinder of radius a = 1 cm deformed into a helix with pitch 2ti /q = 12 cm and helical radius 
h = 0.5 cm rotating with a period T = 1.5 s, the swimming velocity of a helix in a Newtonian 
fluid is U = uqh'^ ~ 0.5 cm/s, which corresponds to unity on the vertical axis of Fig. [3l Note 
that these values are chosen to make aq ^ 0.5, bq ^ 0.25, and h/a ^ 0.5. The time constant 
needed to place the swimming velocity around the maximum value would be on the order of 
0.1 s, which is within the range of time constants found in Boger fluids 26|. Also assuming 
that 1 — ris/i] = rjp/rjs ^ 0.25, again typical for Boger fluids 26|, l3l|], the swimming velocity 
for our example time- varying reciprocal motion is 0. 1-0.5 mm/s. Note that larger swimming 
speeds might be obtained by using fluids with smaller time constant A, but care must be 
taken to remain in the low Reynolds number regime. For the parameters we used above, the 
Reynolds number is less than 0.01 for a Boger fluid with viscosity in the range 0.2-0.5 Pa-s. 



V. CONCLUSION 



In this paper we have calculated the swimming velocity for an infinite cylinder under- 
going arbitrary periodic beating motion in a nonlinearly viscoelastic Oldroyd-B fluid. We 
have assumed that the amplitude of deflection is small compared to the wavelength of de- 
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formations and compared to the radius of the cyhnder. We find that in general, nonhnear 
viscoelastic corrections decrease the swimming velocity relative to the Newtonian case. It 
is interesting to note that the decrease in swimming velocity depends on the nonlinearity: 
in a linear Maxwell fluid, even though the time-reversal invariance of the constitutive rela- 
tion is manifestly broken by memory effects through the time derivative of the stress, for 
a given swimming motion, the flows are exactly the same as the Newtonian case (although 
the stresses are modified) 18|]. Nonlinearity in the constitutive relation is required to rectify 
the changes in the stress into changes of flows. We have used our calculation to examine the 
breakdown of the scallop theorem in a nonlinearly viscoelastic fluid. Within our assump- 
tions of small amplitude deflections of an infinitely long filament, time-reversal- invariant 
motions lead to no net translation. We further showed that net motion in a nonlinearly 
viscoelastic fluid is possible for a more general reciprocal motion, in which the backwards 
and forwards strokes occur at different rates. These predictions could be readily tested 
with simple table-top experiments. The nonlinearity of our constitutive relation necessarily 
limited our analytical results to the case of small deflection. It would be of great value to 
develop a proper local resistance theory for slender rods in a nonlinearly viscoelastic fluid. 

In this paper we have worked in the frequency domain, building arbitrary motions through 
superposition. This allows us to easily treat time-periodic motions, but does not lend itself to 
a physical understanding of motions that suddenly stop or start, which require a complicated 
superposition of many frequencies. In numerical studies of peristaltic flows, Teran, Fauci, 
and Shelley have explored the build-up of stresses from an initial fluid stress given by an 
isotropic pressure. In this context, time-periodic flows correspond to asymptotic long-time 
behavior, which may take a long time to develop in numerical simulations. The numerical 
studies provide a complementary viewpoint well-suited to revealing transient behavior in 
flows. 

This work was supported in part by National Science Foundation Grants Nos. NIRT- 
0404031 (TRP) and DMS-0615919 (TRP), and NIH grant No. ROl GM072004 (CWW). 
TRP thanks the Aspen Center for Physics, where some of this work was done. We thank 
Allen Bower and Eric Lauga for helpful discussions. 
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Appendix A: Forces on a cylinder with prescribed beating pattern. 



The force per unit length acting on the cyhnder is found by integrating the stress tensor 
around the circumference, 

f (s) = a J n(^, z) ■ cr(rb, 0b, Zb)dip, (Al) 

where the dependence of the cyhndrical coordinates {rt, 4>h, 2;b} on ip and s depends on 
the deformation. First consider the force per unit length for a transverse wave with 
polarization in the x-direction only, Eq. Mirror symmetry about the xz-plane implies 
that f±y = 0. To lowest order, fi = r, rb = a, 0b = i^, and = s [see Eq. Q]. Upon 
integration over ip, the force in the z-direction arising from ai^z = Trl^ vanishes, and the 
force per unit length arising from ail^ = — p^^^ + Trr^ and = rj:^ is 

= Rei-^\^ Uhe^^^-^A ] X (A2) 



log a 1 — iDe 

This force is accurate to lowest order in l/log(a). 
Now consider longitudinal traveling waves: 

r{s,'ijj) = ±acostlj + yasim/j + zRe[he''''^'~'^*^ + s]. (A3) 

Using cylindrical coordinates, the position of the surface can also be described by (rb, 0b, -^b)- 
To first order in h, 

rh = a 
0b = ^ 

^b = s + h (A4) 

As in the transverse case, the no-slip boundary conditions at the filament surface are found 
by expanding the argument of the fiow velocity at the surface to first order, 

[v + hz- Vv]„_^_,^, = hz. (A5) 

The nondimensional first and second order fiow equations are the same as in the transverse 
case, but the boundary condition is different. At first order, Vz{a, ip) = —ihz. The first order 
velocity field and pressure are shown in Table [Tll Using these solutions in general expression 
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p 

Vz 

aA 
B 



hKo{r) 
h[aArKo{r)/2 + BKiir)] 


-\h [{B - aA)KQ{r) + aaArKi{r)/2] 
-2Ki{a)/ {a [K^{af - Ki{af] + 2K^{a)Ki{a)] 
aKoia)/ {a [Ko{a)^ - K^{a)^] + 2Ko{a)Kiia)} 



TABLE II: Pressure and velocity field for a cylinder with a longitudinal traveling wave; a = 
(l-iDe)/(l-iDe2). 

for the force per unit length, Eq. (lAip . we find that the first order force for a longitudinal 
traveling wave is 



f(i) 



Re 



2n 1 - iDe2 



-i/ie'(^-*) 



(A6) 



log a 1 — iDe 

As explained before, the first order force is sufficient for calculating the lowest order 
changes to swimming properties that result from shape changes of beating patterns in 
viscoelastic media. To apply Eq. f lA2p . we interpret it in a manner consistent with re- 
sistive force theory [17]. Since the quantity —ihexp{is — it) is the transverse velocity of 
the cylinder, the force per unit length has a complex transverse resistance coefficient of 
— (47r/loga)(l — iDe2)/(l — iDe) for motion occurring with dimensional space- and time- 
dependence exp(igs— iwt). This resistance coefficient depends on q through the dimensionless 
parameter a. As is commonly done in resistive force theory, to treat general motion involving 
a superposition of traveling waves with multiple wavelengths, we replace the slow logarithmic 
dependence on q with a constant involving a typical wavevector (g). The error introduced 
in this process is of order {q — q)/{qlogq). After this replacement we obtain a resistance 
coefficient, in dimensional form, of 

1 - iDe2 



^Oldroyd 



-c 



" 1 - iDe 



(A7) 
(A8) 



log(g)a' 

The parallel resistance coefficient is obtained by a similar analysis for the longitudinal wave, 
and is given by = (^oidroyd^2 g^^^ _ (^^/2. Generalization to motion with multiple 
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frequencies or helical waves may be accomplished by linear superposition. This choice of 
(± and C|| is consistent with those previously introduced for the Newtonian resistive force 
theory. 

With this choice of resistance coefficients, the first-order force per unit length obeys 

f + Af = fvis + ^Afvi„ (A9) 

where fvis is the force per unit length from Newtonian resistive force theory, 

fvis(s,t) = avx(s,t) + C||V||(s,t). (AlO) 

and v_|_(s,t) and V||(s,t) are the perpendicular and parallel components of the fiow velocity 
relative to the rod at the point s on the rod. The velocity v must be interpreted as a relative 
velocity to maintain Galilean invariance. In the limit of zero solvent viscosity, rj^ = 0, these 
equations are equivalent to the resistive force theory proposed by Fulford, Katz, and Powell 
for the linear Maxwell model jisl . 

We can now see why the resistive force theory (IA9I) leads to no change in swimming 
speed relative to Newtonian resistive force theory. First note that in resistive force theory, 
the swimming velocity of a body is determined by demanding that the total force F = / ds f 
on the body vanishes. Integrating Eq. (1A9I) over arclength s and then averaging over a period 
leads to the condition J ds fvis = 0, the same condition for determining the swimming velocity 
in the Newtonian case. 

Although Eq. (1A9|) gives internally consistent predictions for the swimming speed of a 
filament with a traveling wave, it is nevertheless fiawed. To see why, recall our motivation 
of the nonlinear constitutive relation for a fiuid with fading memory in Sec. [Ill There 
we saw that without the nonlinear terms of the upper-convective derivative, the Maxwell 
constitutive relation for an element is sensitive to the local rate of rotation of that element. 
To illustrate this effect, consider shear flow in a channel on a slowly rotating table. A short 
calculation shows that the zero-frequency shear viscosity depends on the rate of rotation of 
the table, an unphysical result if the rate of rotation is low enough that the flctitious forces 



acting on the polymers in solution are small compared to the viscous forces [20|]. 

Likewise, we may consider a slender fllament on a table rotating with speed Q, and ask 
if Eq. (]A9|) predicts forces that are independent of Q. Set ?7s = for convenience. Both the 
fluid and the fllament are rotating with the same speed. Suppose further that the fllament is 
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moved with steady velocity v relative to the fluid. Unless Q is enormously large, we expect 
the force per unit length f to eventually relax to the viscous force per unit length (lAlOl) . 
However, writing Eq. (1A9P in terms the perpendicular and parallel components, and solving 
for the force in the steady state yields 

Thus, we see explicitly that the resistive force theory for the Maxwell constitutive relation 
predicts spurious dependence on the rotation speed fl. Equation (1A9P is only truly reliable 
when working to linear order in deflection. 
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